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Abstract 

The work fluctuations of a harmonic oscillator in contact with a thermostat and driven out of equilibrium by 
an external force are studied experimentally. For the work both the transient and stationary state fluctuation 
theorems hold. The finite time corrections are very different from those of a first order Langevin equation. The 
heat and work fluctuations are studied when a periodic forcing is applied to the oscillator. The importance of the 
choice of the "good work" to compute the free energy from the Jarzinsky equality is discussed. 

Resume 

Etude experimental des fluctuations du travail dans un oscillateur harmonique. Nous etudions 
experimentalement les fluctuations du travail d'une force externe sur un oscillateur harmonique en contact avec 
un bain thermique. Nous trouvons que les theoremes de fluctuations pour les etats transitoires et stationnaires sont 
verifies. Toutefois les corrections de temps fini pour l'etat stationnaire sont tres differentes de celles calculees dans 
le cadre de l'equation de Langevin du premier ordre. Nous etudions aussi le forgage periodique de l'oscillateur. 
Enfin nous discutons l'importance du choix du "bon travail" pour calculer Penergie fibre en partant de l'egalite 
de Jarzinsky. 
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1. Introduction 

Thermal fluctuations play a very important role in the out of equilibrium dynamics of small systems such those 
that one can find in nanotechnologies and biophysics. Indeed in these systems the amount of injected energy 
is often comparable to that of thermal fluctuations, which may produce unwanted and unexpected behaviors. 
For example very large fluctuations may force the system to move in the opposite direction of the one imposed 
by the external forces. In the same way the heat exchanges with the thermal bath are such that the heat may 
instantaneously flow from a cold source to a hot one. Of course these are not violations of the second principle 
of thermodynamics which fixes only the mean values and it does not make any statement for fluctuations. The 
recent fluctuation theorems FTs [f]-[7] allow one to compute the probability of these negative events for the 
work and for the heat flux. It is interesting to notice that the fluctuations of the work done by the external 
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forces to drive the system between two equilibrium states A and B allows one to compute, in some cases, the 
free energy difference AF between A and B, using either the Jarzinsky equality (JE) [9,8] or the Crooks relation 
(CR)[10]. Indeed the JE and CR take advantage of these work fluctuations and relate AF to the probability 
distribution function (PDF) of the work performed on the system to drive it from A to B along any path 7(either 
reversible or irreversible) in the system parameter space. In this paper we will study the possibility of using 
these theorems on real experimental data. Indeed the proofs of FTs and JE are based on a certain number of 
hypothesis; experimenting on real devices is useful not only to check those hypothesis, but also to check whether 
the predicted effects are observable or remain only a theoretical tool. There are not many experimental tests 
of FTs. Some of them are performed in dynamical systems [11] in which the interpretation of the results is 
very difficult. Other experiments are performed on stochastic systems, one on a Brownian particle in a moving 
optical trap [12] and another on electrical circuits driven out of equilibrium by injecting in it a small current [13]. 
The last two systems are described by first order Langevin equations and the results agree with the predictions 
of ref.[7,6].In this paper we study the out-of-equilibrium fluctuations of a harmonic oscillator whose damping is 
mainly produced by the viscosity of the surrounding fluid, which acts as a thermal bath of temperature T. The test 
using a harmonic oscillator is particularly important because the harmonic oscillator is the basis of many physical 
processes [14,15]. The paper is organized as follow. In the next section we will briefly describe the experimental 
set-up. The fluctuation theorem will be studied is section three. In section 4 we will describe the results for the 
work and heat for a the stationary state fluctuation theorem when the oscillator is submitted to a periodic driving. 
In section 5 the JE and Crooks are used to compute the free energy. The importance of the choice of the "good 
work" for the JE is pointed out too. We show indeed that a definition of work which satisfies FT cannot be used 
to compute the free energy with JE. Finally we conclude in section 6. 



2. Experimental set-up 

We recall here only the main features of the experimental set-up, more details can be found in ref.[16,18]. The 
oscillator is a torsion pendulum composed by a brass wire and a glass mirror glued in the middle of this wire. It 
is enclosed in a cell filled by a water-glycerol solution at 60% concentration. The motion of this pendulum can be 
described by a second order Langevin equation: 

he 5 + v + Cd = M + V2k B T v V , (1) 
at* at 

where 6 is the angular displacement of the pendulum, / e ff is the total moment of inertia of the displaced masses, 
v is the oscillator viscous damping, C is the elastic torsional stiffness of the wire, M is the external torque, ks 
the Boltzmann constant and r\ the noise, delta-correlated in time. In our system the measured parameters are 
the stiffness C = 4.5 x 10 -4 Nmrad -1 , the resonant frequency f = \/ C '/ I e g / (2w) = 217Hz and the relaxation 
time r a — 2Ies/v = 9.5ms. The external torque M is applied by means of a tiny electric current J flowing in a 
coil glued behind the mirror. The coil is inside a static magnetic field, hence M oc J. The measurement of 9 is 
performed by a differential interferometer, which uses two laser beams impinging on the pendulum mirror [16,18]. 
The measurement noise is two orders of magnitude smaller than the thermal fluctuations of the pendulum. 6(t) 
is acquired with a resolution of 24 bits at a sampling rate of 8192Hz, which is about 40 times f . The calibration 
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Figure 2. TFT. a) P(W-r) for TFT for various r/r„: 0.31 (o), 1.015 (□), 2.09 (o) and 4.97 (x). Continuous lines arc Gaussian fits, 
b) TFT; S{W T ) computed with the PDF of a). The straight continuous line has slope 1, i.e., S(r) — 1, Vr. 

accuracy of the apparatus, tested at M — using the the Fluctuation Dissipation Theorem, is better than 
3%(see [18]). 



3. Fluctuation theorems 



It is important to make a distinction between Stationary State Fluctuation Theorem (SSFT) and and Transient 
Fluctuation Theorem (TFT) (as defined in ref.[7,6]). TFT implies that the initial state when the force is applied 
to the system is an equilibrium one. In contrast SSFT concerns out of equilibrium stationary states where all 
transients have already relaxed. To study SSFT and TFT we apply to the oscillator a time dependent torque 
M(t) as depicted in Fig. la, and we consider the work W T done by M(t) over a time r: 

ti+T 

Wr = ~ktr I [M{t) ~ M{u)] t dt ■ (2) 

By definition, the TFT implies that U = in eq.2 whereas ti > 3r a for SSFT. As a second choice for M(t), 
the linear ramp with a rising time r r is replaced by a sinusoidal forcing; this leads to a new form of stationary 
state which has never been considered in the context of FT. We examine first the linear forcing M(t) = M t/r r 
(Fig. la)), with M = 10.4 pN.m and r r = 0.1 s = 10.7 r Q . The response of the oscillator to this excitation is 
comparable to the thermal noise amplitude, as can be seen in Fig. lb) where 6{t) is plotted during the same time 
interval of Fig. la). Because of thermal noise the power W T injected into the system (eq.2) is itself a strongly 
fluctuating quantity. 

The FTs state [7] that the probability density functions (PDF) P{W T ) of W T satisfies: 



In 



P(W T ) 



= E(t) W t (3) 



_P(-W T 
where for TFT 

E(r) = 1, Vr (4) 
and for SSFT 

E(t) 1 for r -> oo, (5) 
(see refs.[6,7] for a clear discussion on the difference between TFT and SSFT). 
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Figure 3. SSFT with a ramp forcing, a) PDF of W T for various r/r a : 0.019 (o), 0.31 (□), 2.09 (o) and 4.97 (x). b) Corresponding 
functions S(W T ). c) The slope S(r) of S(W T ) is plotted versus r (□: experimental values; continuous line: theoretical prediction 
from refs. [15,19]) with no adjustable parameter). 



3.1. Transient fluctuation theorem 



We consider first the TFT. The probability density functions P(W T ) of W T are plotted in Fig. 2a) for different 
values of r. We see that the PDFs are Gaussian for all r and the mean value of W T is a few fcgT. We also notice 
that the probability of having negative values of W is rather high for the small r. The function 



S(W T ) = hi 



P(W r ) 



P(-W T ) 



(6) 



is plotted in fig. 2b). It is a linear function of W T for any r, that is S(W T ) = E(r) W T . Within experimental error, 
we measure the slope S(r) = 1. Thus for our harmonic oscillator the TFT is verified for any time r. This was 
expected [2,6], and a derivation of this generic result for a second order Langevin dynamics is given in ref.[15,19]. 



3.2. Stationary state fluctuation theorem 

We now consider the SSFT with tj > 3r a in eq.2. We find that the PDFs of W T , plotted in Fig. 3a), are Gaussian 
with many negative values of W T for short r. The function S(W T ), plotted in Fig. 3b), is still a linear function of 
W r , but, in contrast to TFT, the slope E(r) depends on r. In Fig. 3c) the measured values of E(t) are plotted 
as a function of r. The function E(r) — > 1 for r ^> r Q . Thus SSFT is verified only for large r. The finite time 
corrections of SSFT, which present oscillations whose frequency is close to /o, agree quite well with the theoretical 
prediction computed for a second order Langevin dynamics [15]. We stress that the finite time correction is in this 
case very different from that computed in ref. [6,7] for the first order Langevin equation. The results of Figs. 2,3 
have been checked for several M /t t without noticing any difference. The errors bars in the figures are within the 
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size of the symbols, and they come only from the calibration errors of the harmonic oscillator parameters, and 
statistics of realizations (typically 5 x 10 5 cycles have been used). 



4. Periodic forcing 

We describe in this section the results of the periodic forcing. In this case M(t) = M sin u>dt and the work 
expression (eq.2) is replaced by 



ti+T 



(7) 



where r is a multiple of the period of the driving (r = 2n-K/ui d with n integer). The starting phase ULu d is averaged 
over all possible U to increase statistics. This periodic forcing is a stationary state that has never been studied in 
the details in the context of FT. As already done for the first order Langevin equation [6, 7, 17], we want to consider 
here also the heat fluctuations, that is the energy dissipated by the oscillator towards the heat bath. Multiplying 
cq. (1) by 4| and then integrating between ti and ti + r, we obtain exactly the first law of thermodynamics. The 
change in the internal energy AU T over a time r is: 



AU T = U(U + r) - U(U) = Qr + W T 



(8) 



where Q T is the heat given to the system. Equivalently, (— Q T ) is the heat dissipated by the system. For a harmonic 
oscillator described by a second order Langevin equation (eq.l) the internal energy U(t) and the heat Q T have 
the following expressions: 



U(t) = 



k B T 



-,I e s 



d0(t) 



dt 



(9) 



ti+T 



q t = au t -w t = ~-j^pf I » f t (t') # + j±rf v(t')f t (t'W. 



(10) 



The first term in eq.10 corresponds to the viscous dissipation and is always positive, whereas the second term can 
be interpreted as the work of the thermal noise which have a fluctuating sign. 

We rescale the work W T (the heat Q T ) by the average work (W T ) (the average heat (Q T )) and define: w T = 
T^T ^ T ~ (Q T ) ^ vera S es are time-averages, and they are proportional to r on the stationary state under 
consideration. 

As an example we consider a measurement performed at M — 0.78pN m and at ajd/(27r) = 64Hz. The PDFs of 
w T , AU t and q T are plotted in Fig. 4 for different values of n. The average of AU T is clearly zero because the time 
r is a multiple of the period of the forcing. The PDFs of the work (fig. 4a) are Gaussian for any n whereas the 
PDFs of heat fluctuations q T have exponential tails(fig. 4c). These exponential PDFs can be understood taking 
into account that, from eq. (10), — Q T = W T — AU T and that AU T has an exponential PDF independant of n 
(fig. 4b). In the case of the periodic forcing the PDF of q T can be exactly computed (see ref.[19]): 



f (Qr) 



■ exp 



2 



,Qt~{Qt) 



erfc 



- (Qt) + Ow 



+ e 



-(Q t -(Qt» 



erfc 



}t + (Qr) + Ow 



,(H) 



where erfc stands for the complementary Erf function and o w is the work variance. In Fig. 4c, we have plotted 
the analytical PDF from eq. (11) together with the experimental ones, using the values of Ow and (Q T ) from the 
experiment. The agreement is perfect for all values of n and with no adjustable parameter, using eq. (11). 
To quantify the symmetry of the PDF around the origin, we define the function S as: 



S(e T ) ee 



(Et) 



In 



' P{e T ) - 



(12) 
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Figure 4. Sinusoidal forcing, a) PDFs of w T . b) PDFs of AUt- c) PDFs of q T . d) Functions S(q-r). In all plots, the integration time 
r is a multiple of the period of forcing, t — In-njuJd-, with n — 7 (o), n — 15 (□), n — 25 (o) and n — 50 (x). Continuous lines 
are theoretical predictions with no adjustable parameter, c) The slope S q (n) of S(q r ) for q T < 1, plotted as a function of n (o). 
The slope Stu(n) of S(io T ) plotted as a function of n(d). Continuous line is theoretical prediction, f) The slopes E(n), plotted as a 
function of n for two different driving frequencies uJd — 64 Hz (□) and 256 Hz (o); continuous lines are theoretical predictions from 
ref.[15,19] with no adjustable parameter. 
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where e T stands for either w T or q T and E T stands for either W T or Q T . The question we ask is whether: 



lim S(e T ) = e T 



(13) 



as required by SSFT. S(q T ) is plotted in Fig. 4d) for different values of n ; three regions appear: 

(i) For large fluctuations q T , S(q T ) equals 2. When r tends to infinity, this region spans from q T — 3 to infinity. 

(ii) For small fluctuations q T , S(q T ) is a linear function of q T . We then define E 9 (n) as the slope of the function 
S(q T ), i.e. S(q T ) = E 9 (n)g T . This slope is plotted in Fig. 4e where we see that it tends to 1 when t is 
increased. So, SSFT holds in this region II which spans from q T = up to q T = 1 for large r. 

(iii) A smooth connection between the two behaviors. 

The PDF of the work being Gaussian, the functions S(w T ) are proportional to w T for any r, i.e. S(w T ) = 
Y, w {n)w T (ref. [15]). T, w (n) is plotted in Fig. 4e) and we observe that it matches experimentally E,(n), for all 
values of n. So the finite time corrections to the FT for the heat are the same than the ones of FT for work [15]: 

E w (n) = E,(n) = 1 + K/n + 1/n O (e~ Tn/Ta \ where K is a constant. 



We find that that the scenario plotted in fig.4a)-e) is the same for any driving frequency ujd [15,19]. However 
the finite time corrections of E m (n) are a function of u)d- The function E w (n), measured at uJd/^ = 64Hz, is 
compared with that measured at tUd/2-K — 256Hz in in fig.4f). The convergence rate is quite different in the two 
cases, which agree with the theoretical predictions for a second order Langevin equation (see ref. [15,19]). 

We also notice that in the case of the sinusoidal forcing, the convergence is very slow in Figs. 4e) and 4f), 
we see that it takes several dozens of excitation period (500 ms for uJd/^ = 64Hz) to get E(n) = 1 within one 
percent. On the contrary for a ramp forcing, this is achieved in about (20 ms), i.e. a few r a (see Fig. 3c). 

In this section we have seen that also in the case of the sinusoidal forcing the agreement between the computed 
and measured finite time corrections is very good. These results prove not only that FTs asymptotically hold for 
any kind of forcing, but also that finite time corrections strongly depend on the specific dynamics. The detailed 
theoretical analysis of this behavior for a second order Langevin equation can be found in ref. [19]. 



5. TFT and the Jarzynski equality the Crooks relation 

In this section we want to discuss the relationship between TFT and JE. The Jarzynski equality (JE) relates the 
free energy difference AF between two equilibrium states A and B of a system in contact with a heat reservoir to 
the PDF of the work performed on the system to drive it from A to B along any path 7 in the system parameter 
space. Actually thermal fluctuations are here very useful because they allow the system to explore all possible 
paths from A to B. Specifically, let us consider the work done by the external torque on the harmonic oscillator 
described in section 3. The external torque M(t) drives the system from M(0) = (equilibrium state A) to 
M(t > T r ) — M (equilibrium state B). Jarzynski defines the work performed by M(t) to drive the system from 
A to B as 



is the classical work such that — W cl / {UbT) is equal to W Tr defined in eq.2. Here we define W° such that it has 
the same sign of W J . Thus W J and W cl are related but they are not the same and they differ by a boundary 
therm, which makes an important difference in the fluctuations of these two quantities. This difference between 
the W and W c i, has been already pointed out in ref. [20]. One can consider an ensemble of realizations of the 
"switching process" from A to B with initial conditions all starting in the same initial equilibrium state. Then the 
work may be computed for each trajectory in the ensemble. The JE states that [9] 




(14) 



where 




(15) 







AF = -\ log(exp [-(3W J \), 



(16) 
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Figure 5. Mo = 11.9pNm (a) P f (W) and P b (-W), (b) P f (ff cl ) and P b (-W cl ); M = 6.1j>JVm: (c) P t (W) and P b (-W), (d) 
Pf(W cl ) and Pb( — W^ 01 ) (experimental forward and backward pdfs are represented by o and □ respectively, whereas the continuous 
lines are Gaussian fits) 

where (•} denotes the ensemble average, (i~ x = ksT. In other words (exp [— /JWdiss]) = 1, since we can always 
write W = AF + Wdiss where Wdiss is the dissipated work. In our experiment we can also check the CR which is 
related to the JE and which gives useful and complementary information on the dissipated work. Crooks considers 
the forward work Wt to drive the system from A to B and the backward work W b to drive it from B to A. If the 
work pdfs during the forward and backward processes are Pf(W J ) and Pt>(— W J ), one has [10,8] 



Pt(W J ) 
P b (-W J ) 



exp (f3[W J - AF]) = exp [f3W diss ]. (17) 



A simple calculation from Eq. (17) leads to Eq. (16). However, from an experimental point of view this relation 
is extremely useful because one immediately sees that the crossing point of the two pdfs, that is the point where 
P f (W J ) = P b (-W J ), is precisely AF. Thus one has another mean to check the computed free energy by looking 
at the pdfs crossing point W x ■ 

Many details of the test of the JE and CR in the harmonic oscillator can be found in ref. [18]. Here we want 
only to stress the difference between the fluctuations of the classical and the Jarzinsky works. This can be seen 
in fig. 5 where we plot the the pdfs of W J (figs. 5 a,c) and W cl (figs. 5 b,d) for the forward and backward processes 
measured for two values of M specifically M — ll.QpNm (figs.5a,b) and M = G.lpNm (figs.5a,b). The stiffness 
C of the oscillator was 7.5 10~ 4 Nm rad^ 1 in this experiment. In Figs. 5, bullets are the experimental data and 
the continuous lines their fitted Gaussian pdfs. Comparing the pdfs of W J and W cl for the same M a we see 
that the variance of W cl is much larger than that of W J . Indeed W cl presents both positive and negative values 
and, as we have seen in sect. 3, P(—W cl ) satisfies TFT. Thus it is clear that if W J satisfies eq.16 this equation 
cannot be satisfied by W cl . Indeed using the measured values of W J in eq.16 one finds [3AF — —23.3 ± 0.4 for 

M = 11.9 and /3AF — —6.3 ± 0.3 for M = 6.1. These values correspond to AF = —-^q for the driven harmonic 
oscillator. In Fig. 5, the pdfs P { (W) and P b (-W) cross at (3W ~ -23.5 (fig.5a) and in /3W ~ -6.1 (5c), which 



again correspond to the expected values. Thus we see that JE and CR can be safely applied on experimental 
data to measure AF if the work W J defined in eq.14 is used. In ref.[18] it has been shown that this result is true 
independently of the ratio T r /r a and of the maximum amplitude of M a . This has been checked at the largest M a 
and the shortest rising time r r allowed by our apparatus. 

The fact that the JE and CR are satisfied only by W J but not by W cl can be easily understood by noticing 
that a variable that satisfies TFT cannot be used in eq.16. Indeed one may rewrite eq.16 as: 



oc 

cxp-/3AF = < exp(-/3 W) >= j exp(-f3 W)P(W)dW (18) 

— oo 

where W is an energy injected into the system on a time r. If P(W) satisfies TFT then from eq.3 and eq.4 : 

P(W) = exp(-/3 W)P(-W), Vr 
. Replacing this identity in eq.18 one finds: 

oc oc 

J exp(-/3 W) P(W) dW = j exp(-/3 W) P{-W) exp(/3W) dW = 1 (19) 

— oc — oo 

showing that if W satisfies TFT then cannot be used to compute AF. Thus one concludes that as W cl satisfies 
TFT it cannot satisfy the JE as it shown in eq.19. This observation points out the importance of the choice of the 
"good variables" in order to use the JE to compute the free energy difference in a more complex system. Several 
other limitations for a safe application of JE and CR on real data have been already discussed in ref.[18]. 



6. Conclusions 



In conclusion we have applied the FTs and JE to the work fluctuations of an oscillator driven out of equilibrium 
by an external force. The experimental results agree with theoretical predictions both for TFT and SSFT (eqs.3,4 
and 5). Indeed we observe that in the transient case when the system starts in equilibrium (TFT) the FT (eq.3) 
holds for any time, whereas when the system is in a stationary out of equilibrium regime (SSFT), eq.3 holds only 
asymptotically. The SFFT presents a complex convergence to the asymptotic behavior which strongly depends 
on the form of the driving. The exact formula of this convergence can be computed using several experimental 
evidences of the statistics of the fluctuation. The details of derivation of the finite time correction and of the PDF 
for the heat will be the subject of a long paper [19]. We have also discussed the importance of the choice of the 
variable in order to safely apply JE. The results reported in this paper are useful for many applications going 
from biological systems to nanotechnology, where the harmonic oscillator is the simplest building block. 
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